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Abstract 

There are (at least) two surfaces of particular interest in eternal de Sitter space. 
One is the timelike hypcrsurfacc constituting the lab wall of a static patch observer 
and the other is the future boundary of global de Sitter space. We study both linear 
and non-linear deformations of four-dimensional de Sitter space which obey the Ein- 
stein equation. Our deformations leave the induced conformal metric and trace of the 
extrinsic curvature unchanged for a fixed hypersurface. This hypersurface is either 
timelike within the static patch or spacelike in the future diamond. We require the 
deformations to be regular at the future horizon of the static patch observer. For lin- 
earized perturbations in the future diamond, this corresponds to imposing incoming 
flux solely from the future horizon of a single static patch observer. When the slices 
are arbitrarily close to the cosmological horizon, the finite deformations are character- 
ized by solutions to the incompressible Navier-Stokes equation for both spacelike and 
timelike hyper surfaces. We then study, at the level of linearized gravity, the change 
in the discrete dispersion relation as we push the timelike hypersurface toward the 
worldline of the static patch. Finally, we study the spectrum of linearized solutions 
as the spacelike slices are pushed to future infinity and relate our calculations to 
analogous ones in the context of massless topological black holes in AdS4. 
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1 Introduction — X+/Static Patch Schizophrenia 

The observation that the universe is dominated in energy by a positive cosmological 
constant [U [2] suggests that it may evolve into an asymptotically de Sitter era in the 
future. In such a context, any surviving observers will find themselves surrounded by a 
cosmological horizon whose size is set by the value of the cosmological constant. The 
nature of this horizon is rather enigmatic given that an observer can never reach her 
cosmological horizon - it is always at a finite fixed distance away from her. Thus, unlike 
a black hole, the static patch horizon is not localized in some finite region of space. On 
the other hand, it behaves classically, and to an extent quantum mechanically, very much 
like a black hole horizon. For instance there is a temperature and entropy associated with 
the cosmological horizon [3]. 
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One of the main challenges confronted by theorists is to uncover the nature of the 
holographic principle in the context of asymptotically de Sitter universes. One may be 
inclined to propose that de Sitter holography should only describe a single patch, given 
that this is the region of space accessible to a single observer [H O |6l [71 [H |9l [lOl [TTl [121 
[T3t [T^ I15j . In fact, the observer's worldline in the static patch resembles in many ways 
the boundary of anti-de Sitter space (in the presence of an eternal black hole) and may 
constitute the ultraviolet regime of the fundamental description [I2l[l31|15], perhaps in 
a way related to matrix theory. In such a case, the dynamics of the cosmological horizon 
would constitute the deep infrared behavior of the putative worldline theory very much 
like the near horizon dynamics of a black hole in anti-de Sitter space constitutes the deep 
infrared behavior of the dual conformal field theory. Such behavior rather universally 
takes the form of fluid dynamics. 

There have been several attempts to relate general relativity to fluid mechanics dat- 
ing back to the 1970s with the black hole membrane paradigm \T7\ [T8] (see [19] for 
an application to de Sitter space). The membrane paradigm focuses on the observation 
that the equations governing the dynamics of horizon surfaces in general relativity can 
be written in a form analogous to that of the Navier-Stokes equation of fluid mechan- 
ics. However, whilst finding a striking analogy, the central equation of the membrane 
paradigm is often referred to as the Damour-Navier-Stokes equation, highlighting the 
fact that it differs from the Navier-Stokes equation in key ways. Building on this, recent 
papers |20| [2T| [22] constructed a setup where near horizon dynamics in gravity precisely 
relates the Einstein equation to the incompressible Navier-Stokes equation. These stud- 
ies were also inspired by analyses of connections between gravity and fluid mechanics 
in the context of the AdS / CFT correspondence |23[ [2^ [23] and the low energy limit of 
the dual field theory. Given the striking similarities between the thermodynamics of a 
black hole horizon and a cosmological horizon, it is natural to extend such a fluid/gravity 
correspondence to include spacetimes with a cosmological horizon. 

Yet another natural location for a definition of quantum gravity in de Sitter space is 
the spacelike boundary at future infinity, known as X+ [261 [271 [28l [291 [301 [HJ [32] . The 
'metaobservables' on are given by correlators between causally disconnected points. 
In a sense, we are metaobservers of the 'would be' X"*" of the inflationary de Sitter era. 
Such metaobservables are conjectured to constitute the correlation functions of a lower 
dimensional non-unitary Euclidean conformal field theory. Recently, an exact example of 
this dS/CFT correspondence was proposed in [33]. It is not at all clear how the physics of 
the static patch observer is captured by the theory at X+. The static patch observer can 
at most observe a single point (or tiny region) of where her worldline intersects the 
future boundary. Thus, from the static patch observer's point of view, fixing the geometry 
outside her future horizon is akin to fixing a gauge since it will never affect the physics 
she observes |14j . In particular, one might envision fixing the geometry near by fixing 
the data on a spacelike slice in the future diamond and allowing only flux originating 
from a single static patch to come through. This is somewhat analogous to the boundary 
condition that there is no incoming flux from the past horizon of a black hole. It may then 
be speculated that the existence of a finite number of 'holographic reconstructions' of a 
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static patch observer in the dual CFT at X"*" would be a manifestation of the finiteness 
of the de Sitter entropyj^ 

After reviewing the classical geometry of de Sitter space, the first part of this paper 
will explore some of the classical features of the cosmological horizon as viewed by an 
observer in a purely de Sitter universe - the static patch observer. We examine the 
Einstein equation both linearly and non-linearly and uncover that the solutions are char- 
acterized by solutions to the incompressible Navier-Stokes equation on a two-spherej^ 
This same equation recently appeared in the context of the Schwarzschild black hole [22j 
and requires the velocity field Vi{T,n^) where ri* = {9,(j)} and the pressure P{t,0,^) to 
satisfy 

drv' + V's2P + VjVl.,v' - V (V\,v' + Ry) = , V^^Vi = (1.1) 

where v is the viscosity. Indices are raised and lowered with respect to the round metric gij 
on the S"^ of radius rs for which Rij{= gij/f'g) is the Ricci tensor. At the linearized level, 
this is done by imposing Dirichlet boundary conditions on a timelike surface arbitrarily 
close to the cosmological horizon and the absence of incoming flux from the past horizon of 
the static patch. These boundary conditions resemble the solipsistic boundary conditions 
of |15j . which allow for an examination of the isolated static patch dynamics, unperturbed 
by external sources from the past horizon. We find that the linearized solutions must 
obey the dispersion relation of the incompressible, linearized (pressureless) Navier-Stokes 



equation (3.13). At the non- linear level, again in a near cosmological horizon expansion, 
we impose (conformal) Dirichlet boundary conditions on a timelike slice and regularity 
of the solutions as they approach the future horizon. By (conformal) Dirichlet boundary 
conditions, we mean analysing perturbations which leave the induced geometry on a 
fixed timelike hypersurface of constant extrinsic curvature unchanged up to a conformal 
factor]^ Then, we comment briefly on the possibilities of deforming this non-linear fluid 
by placing a small black hole at the origin of the static patch. In an attempt to connect 
our fluid dynamical modes to the analogous excitations of the worldline, which are the 
quasinormal modes, we return to the linearized analysis to study how the linearized 
dispersion relation varies as we push the surface from the cosmological horizon to the 
worldline. 

In the second part of the paper we make some mathematical observations about space- 
like slices foliating the region outside the future horizon of the static patch. We examine 
the behavior of linearized solutions to the Einstein equation near, but outside, the future 
cosmological horizon. Our solutions are subjected to Dirichlet boundary conditions on 
a fixed spacelike surface and to contain incoming flux solely from a single static patch 
observer. We find a discrete set of modes obeying the dispersion relation of the linearized 
Navier-Stokes equation, where the time coordinate has become the non-compact spacelike 



^This may ultimately be related to the mysterious use of the Cardy formula for 'counting' the de Sitter 
entropy in some lower dimensional examples |34l I35| . 

■^As in we analyze the metric through the first three orders in a near-horizon expansion. A 

generalization of the all-orders proof of 36 might be possible in our case, but we will not attempt to do 
so herein. 

^Henceforth, in the non-linear analysis, we will refer to these boundary conditions as Dirichlet bound- 
ary conditions. 
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coordinate moving us along the spacelike slice. The non-linear solutions to the Einstein 
equation which satisfy Dirichlet boundary conditions on the spacelike slice and which are 
regular at the horizon from which flux is coming, are indeed characterized by solutions to 
the incompressible Navier-Stokes equation. The Navier-Stokes equation uncovered here 
on the spacelike slice is equivalent to that discussed in the context of the timelike surface, 
except that the sign of the viscosity is flipped. We end by noting that the setup of the 
problem in this future diamond of de Sitter space, and in particular the pole structure at 
, is connected by an analytic continuation to analogous problems in Lorentzian AdS4 
with hyperbolic slicing. 

2 Geometry and Framework 

In what follows we will study the geometries of several patches of de Sitter space pertinent 
to our analysis. Instead of the global patch of de Sitter space containing the past and 
future infinities, denoted by X~ and X"*", we will focus on patches that are more suited to 
the description of local observers. 

2.1 The Static Patch 

The four-dimensional static patch metric solves the Einstein equation in the presence of 
a cosmological constant A > 0, 

= G/.^ + A5m^ = (2.1) 

and is given by: 

ds^ = - (1 - {r/if) df + (1 - {r/lfy^ dr'^ + r^dn^ , (2.2) 

where r G [0, £], t G M and is the round metric on 5^. The quantity i is the de Sitter 
length and is related to the cosmological constant as A = -|-3/£^. The above metric covers 
a quarter of the global de Sitter geometry, it describes the intersection of the future and 
past causal diamonds of a constant r worldline beginning at I~ and ending at We 
call this the Southern patch of de Sitter space. 

One notices that r = i corresponds to a cosmological event horizon, beyond which 
events are forever out of causal contact from the Southern observer. The Killing vector 
dt becomes null at r = i and the above coordinate system breaks down. 

The Southern patch can be smoothly connected to another region covering an addi- 
tional quarter of de Sitter space, by continuing the above metric to r G [i, oo]. For r > i, 
t becomes a spacelike coordinate and r becomes timelike. We can consider gluing two 
such regions, one behind the past cosmological horizon, known as the past diamond con- 
taining I~ , and the other beyond the future cosmological horizon, known as the future 
diamond containing X"*" . 

The remaining quarter of the global de Sitter space is given by an additional static 
patch system known as the Northern patch. The Southern and Northern patches each 
intersects at a single point. In figure [l] we demonstrate the several patches discussed 
above in a Penrose diagram. 
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Fig. 1: Penrose diagram of de Sitter space indicating the various static patches and 
future/past diamonds. 



2.2 Null Foliations 

It will be convenient to introduce an additional coordinate system which smoothly covers 
both the Southern patch and the future diamond. This is achieved by the following 
coordinate transformation: 

dv T 

= - (1 - {r/m ' " = ^'-^^ 

leading to the metric 

= -{l-v^)du'^ -2dudv + v'^d^l . (2.4) 

Up to a constant time shift we find u£ = t — £tanh~^ r/£. Constant u surfaces are null 
lines emanating from the origin at f = and ending at X"*" where v = oo. The norm of 
the Killing vector du changes sign at v = 1. 



2.3 Approaching the Horizon 

Finally, we would like to introduce a dimensionless parameter a > allowing us to 
approach the cosmological horizon. In order to achieve this, we rescale time to u = T/2a 
and define p = {1 — v)/2a. As we take the limit a — )• 0, we redshift time and for any 
finite p, v will be forced to lie near the cosmological horizon. The metric is given by: 

^ = (-P+pA dr"^ + 2dTdp + (1 - 2apf d^l . (2.5) 

The coordinate range of p covering the Southern patch is given by p G [0, l/2a] and the 
norm of du vanishes at p = 0. The constant p and r surfaces are shown in figure [2} 



6 




Fig. 2: Penrose diagram of de Sitter space indicating constant p (red) and r (gray diag- 
onal) slices. 

As opposed to the Scliwarzschild case, where a similar expansion would continue for 
indefinite powers of a, the above expansion terminates at order 0{q^). This is due to 
the absence of a term ~ 2M/r in the Qrr component. We could of course add such a 
term, which would correspond to introducing a small mass or black hole centered at the 
origin of the static patch. 

3 Incompressible Fluids 

Having specified the geometry relevant to our problem, we proceed to discuss the nature 
of perturbations solving the Einstein equation with positive A near the cosmological 
horizon. We begin with a linearized analysis. 

3.1 Linearized Analysis 

Linearized gravity about spherically symmetric spaces with non-zero cosmological con- 
stant was examined in |371 138j . The two gravitational degrees of freedom transform as 
a (divergenceless) vector and a scalar under the SO{?>) symmetry of the S"^ . There is 
no transverse-traceless tensorial spherical harmonic for a two-sphere. Let us work in a 
gauge where 5gij = for G {^}- The metric vector perturbations can be expressed as: 

5git = X (1 - (r/^)2) (1 + rdr) , (3.1) 
^9ir = X (T^^Jt^^*'^- (3-2) 
The vector spherical harmonic Vi satisfies the following relations on the unit two-sphere: 

(V|2 + fc^) Vi = , V^aVi = , (3.3) 
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with eigenvalues are ky = /(/ + 1) — 1 and / = 1,2,... The master field $t, obeys the 
master equation: 

'V?„,-'<^V. = 0, (3.4) 



where g^"^^ corresponds to the two-dimensional de Sitter static patch. A similar result 
holds for the scalar perturbations, which we discuss in appendix [K\ 

The solutions to the above equation were analyzed in [Tl] and are found to be hy- 
pergeometric functions. For our purposes we would like to obtain the linearized solu- 
tions in the null coordinate system (2.4). Assuming a Fourier decomposition in time, 

= e~'^^°^'^^^'^'^^^^ipv{p), the equation of motion becomes: 

; (1 - apf dl + 4p (1 - ap) (1 - 2ap) dp+ 

= . (3.5) 



4a2(l - 2apfuj'^ - Aap (1 - ap) {k^ + l) 



(1 - 2apY 

The two linearly independent solutions for / > 1 are given by: 

^oni ^ 2i^i[ai,6i;ci;ap(-l + 2ap)-i] (1 - ap)"*"" (1 - 2a/,f ,(3.6) 

^ p+iau J^^^ ^ 2ap)-i] (1 - ap)-^"" , (3.7) 

with: 

ai = — / — 2iau} , bi = 1 + I — 2iauj , ci = 1 — 2iauj ; (3-8) 

02 = —I , b2 = 1 + I , C2 = 1 + 2iaijj . (3-9) 

The superscripts 'out' and 'in' indicate that the mode is purely outgoing at the future 
horizon or purely incoming from the past horizon. The above expressions are linearly 
independent solutions for (cj — Ui — hi) = 2iaijj ^ Z (see [39]). In the case where (cj — 
o-i — bi) = 2iauj is an integer, logarithmic solutions will appear. Given that a2 and 62 
are integers c/J^" is in fact a finite polynomial for 2iauj ^ Z, as it can be shown that the 
hypergeometric series terminates. For / = 1, the linearized perturbations become time 
independent and are like the introduction of a small amount of angular momentum (we 
discuss this case in appendix [B]) . 

The linearized purely outgoing metric components (3.1) in the (r, p)-coordinate sys- 
tem become: 



X„out r)'^ ; ^^Q'^+l/i „ „\—iauj+l ( -\ 2(y,p) \ nut /o i ^^ 

(>giT =2ViXe p ^ [l-ap) ^ (1 — dp j ip^ , (3.10) 



out 



bg'^' = -2aVi X e' 



1 — ap 
P 



^ _ iujil-2ap) \ _ {l-2ap) g ' 
2p{l-ap) J 2a 



, „out 



(3.11) 

Both and bg°^^ are regular at the future horizon p = 0. 
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3.2 Linearized Fluid Modes 



Having written down the linearized solutions, we now discuss the choice of boundary 
conditions. We impose Dirichlet boundary conditions on a given timelike hypersurface 
at some fixed p, and without loss of generality, we choose the location of this timelike 
hypersurface to be at /? = 1. Taking a — )• pushes this hypersurface arbitrarily near the 
cosmological horizon and thus allows us to probe the near horizon dynamics. 

Our particular Dirichlet boundary condition, shown in figure [3j is that the linearized 
perturbations are purely outgoing and leave the intrinsic geometry of the p = 1 hyper- 
surface unchangedrl Imposing 6g"^''{p = 1) = enforces a discrete dispersion relation, 




Fig. 3: Our boundary conditions for the linearized modes are such that the induced metric 
on the /O = 1 slice is unchanged and there is no incoming flux from the past horizon. 

which to leading order in a is given by: 

ujf = + -2) , 1 = 1,2,... (3.12) 

We interpret these linearized modes as fiuid modes of the velocity field of the incom- 
pressible, linearized (pressureless) Navier-Stokes equation on a sphere: 

drVi = V {Vl2Vi + RijV^) , V's2Vi = , (3.13) 

where the viscosity u = 1. The incompressibility of the fluid is equivalent to the vanishing 
divergence of Vj, which can be seen by identifying Vi ~ e~^^'^Vi. We further note that the 



explicit modes (3.10) with uj = Uf decay in time and are regular at the future horizon 

p = oj3 

By arguments similar to those in |20j . one expects that the result v = 1 corresponds 
to a ratio of shear viscosity to entropy density which is l/47r. This suggests that the 



*This is the simplest choice of Dirichlet boundary conditions and thus allows for a clear analysis. In 
general, we can choose more involved Dirichlet boundary conditions on the p = 1 hypersurface. 

^If these modes are taken back in time to t — >■ — oo they diverge and the perturbative solution is no 
longer reliable. As usual we only consider wavepackets of the linearized solutions which are finite for all 
asymptotia. 
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incompressible fluid we have found near the de Sitter horizon shares this feature with 
the fluids found near the Schwarzschild, Rindler and planar AdS black hole horizons 
[201 ED Ea HO] (see also [E]). 

There is also a decoupled set of scalar excitations which transform as scalars under 
the SO{3) of the S'^. However, the incompressibility condition implies that we could only 
consider the spherically symmetric scalar mode, which reduces the fluid vector fleld to a 
trivial one(see appendix [A|) . We will not consider such modes in what follows and simply 
set them to zero in the linearized analysis. 



+ a 



3.3 Non-linear Analysis 

Having analyzed the linearized case, we now turn to the question of non-linear deforma- 
tions. The analysis follows directly the Schwarzschild case analyzed in [22j|^ In particular, 
in this non-linear analysis, we impose the (conformal) Dirichlet boundary conditions on 
the hypersurface as in |22j . 

To be precise, we will consider the following finite deformation of the static patch 
geometry as an expansion in aj^ 

^ ^ (3.14, 

+ p^dr^ + 2dTdp + dn^ + (1 - p) [u^dr^ - 2v^dTdx'] - IpPdr'^ 
{-ip + 2P) d^ll + (1 - p)v^Vjdx'dx^ 

- [p^ - l) (V\i + RjiV^) drdx' - 2vidpdx' + {v"^ + 2P) drdp + 2(1-/9) tjjfUTdx' 

+ o? (Ap'dnl + 2glf Ux'dp + g^fUx'dx^^ +.... 

The t;*, P and <^f^ are functions of (r, fi*) only while the gl^ and g^^ are functions 
of (r, p, il*). We have chosen a gauge where gpp = 0. As boundary conditions we require 
the perturbations to preserve the induced metric on the hypersurface p = 1 

ds\^ = (^-^ + 1^ dr"^ + (1 - 2a)^ dQ^ , (3.15) 

up to a conformal factor 

l + 2aP + 0{a'^) . (3.16) 

We also study perturbations such that this hypersurface has constant mean extrinsic 
curvature and that the solution is regular at the future horizon p = 0. These boundary 

®It should be noted that we have presented a more complete linearized analysis than would be possible 
for the Schwarzschild case, given t he ex istence of exact linearized solutions in dS4. 

^When writing out the metric (3.141 we have omitted the metric components grr^ffTT ,gpT ,<?}t 
and higher order contributions since these do not affect the Einstein equation to the order that we 
consider. 
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conditions are the natural extension of the boundary conditions imposed on the linearized 
fluid modes of the last section. 

We now examine the conditions on the deformation parameters imposed by the Ein- 
stein equation with a positive cosmological constant Qfj,i, = up to and including 0{a^). 
We further assume that the only excited field is the metric. The first non-trivial condition 
appears at 0{a~^). Here, for G^-r = to be satisfied, the velocity field is required to 
be incompressible. At the next order 0{a^), the non-trivial equations are Qtt = Gri = 0. 
From the Gri = 0, it follows that the (u*, P) need to satisfy the non- linear incompressible 



Navier-Stokes equation (1.1) on a unit S'^. Our result is in complete accordance with the 
linearized analysis]^ 

From the Gtt = Einstein equation at 0{a^) we find the requirement 

V^a^"^ = 2drP + driv"^) + total derivatives on the 2-sphere. (3.17) 
The above relation implies: 

dr {^j v^dn2^ = -2dr {^j PdVt^ . (3.18) 

There is a similar condition between the velocity and pressure fields in the Schwarzschild 
case [22jP| This component of the Einstein equation also determines a scalar function 

mvolvmg gy \g\- ' . 

3.4 Deformations of the Fluid 

A natural question to ask about the fluid is whether one can deform it. In this section, 
we discuss two simple examples of possible deformations of the fluid. 

The first is given by adding a small non-rotating black hole of mass M at the origin 



of the static patch. This changes the —gtt = 9" components of the metric (2.2) to 
V{r) = 1 — {r/l)"^ — 2M/r. For positive values of M, adding the black hole has the effect 
of pulling in the cosmological horizon and thus decreasing its size. For small e = M/t, 
the new position of the cosmological horizon is given by rcos = ^(1 — e) to leading order. 
In the analogous case of the Schwarzschild black hole, placing a mass at the center of the 
static patch corresponds to extracting some mass from the black hole, thus shrinking its 
horizon. The mass deformation we have described preserves the spherical symmetry of 
the background and thus the near horizon dynamics are expected to be governed by the 
Navier-Stokes equation on a sphere. 

A slightly more involved deformation corresponds to placing a small rotating mass at 
the origin of the static patch. This will cause the cosmological horizon itself to rotate. 



We would not expect to see the pressure in the linearized analysis since at the linear level vector 
and scalar representations of 5*0(3) decouple. This is no longer the case at second order in perturbation 
theory where we expect the equation for the vector representation to be affected by scalars as in the non- 
linear case. As in the linearized analysis, we have not considered sound modes, which would contribute 
to the divergence of v''. 

^It has been speculated that such an integral relation is related to changes in the horizon area [221153] . 
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The function determining the positions of the horizons is now given by: 

V{r) = (1 + {a/rf) (l - {r/lf) - 2M/r . (3.19) 

As with the mass term, adding angular momentum shrinks the cosmological horizon. To 
lowest order in small e and small v = {a/t)^ we find that: 

rcos = t{l-e{l-v)-^e^ + 0{£\£'v,ev\v'^)^ . (3.20) 

The angular momentum of the space-time becomes Qq^ = —aMj (l + [ajPf'^. It should 
be noted that a finite deformation with angular momentum will also deform the sphere 
into a spheroidal surface. Thus we lose spherical symmetry and it might be possible that 
the near horizon dynamics is no longer governed by the Navier-Stokes equation on the 
round metric of S"^ . 



4 Pushing the TimeUke Surface 

So far we have analyzed the behavior near the cosmological horizon. Another timelike 
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surface of interest in the static patch is given by the observer's worldline at r = 
Returning to the analysis of linearized gravity, if we impose Dirichlet boundary conditions 
leaving the worldline unperturbed for purely outgoing modes, we obtain another set of 
discrete modes known as quasinormal modes (see for example |41j ) . In the original static 



patch coordinates (2.2) these are given (for the vector modes) by: 

tj„£= -i(n + / + l) , n = 0,1,2,... (4.1) 
Due to the fact that Z > 1, a gapless mode is absent in the above spectrum of quasinormal 



modes. This is in contrast to the fluid modes (3.12) which have = at / = 1. The 
gapless mode is absent due to the fact that the w = 0, I = 1 perturbation diverges on the 
worldline, as shown in appendix [B] It reappears in the spectrum as soon as we 'puff up' 
the thickness of the worldline. 

4.1 'Flowing' the Dispersion Relations 

Our aim is to study the behavior of perturbative data on constant r surfaces as we push 
them from the horizon toward the worldline. There is a clear distinction between the 



lowest n = quasinormal modes (4.1) and the fluid modes (3.12). Given a constant r 
slice at some position r = Tc we impose Dirichlet boundary conditions leaving the induced 
metric on the r = Tc unchanged. This constant r surface is directly analogous to the 
timelike hypersurface at /) = 1 considered above. Furthermore, we require that the modes 



'^"Due to its resemblance with the boundary of AdS in the presence of an eternal black hole, recent work 
has emphasized the potential importance of the worldline as a candidate for the ultraviolet (holographic) 



description of the static patch [12lll5j . Although such a holographic duality is far from clear, one expects 
that the infrared behavior must give rise to the Navier-Stokes equation described in the former section. 
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are purely outgoing. As before, these two conditions will only be satisfied for a discrete set 
of modes, but the dispersion relation will now depend on the dimensionless parameter 
X = Tc/I!.. For the surface near the horizon we have x — )• 1 and as we approach the 
worldline we have x — )• 0. For general x, the problem cannot be approached analytically 
and we must resort to numerics. 

The dispersion relation corresponds to the pole structure of the Green's function of 
the vector modes on the particular cutoff surface. Thus, naturally, a flow of the dispersion 
relation corresponds to a flow of the Green's function itself. For an incompressible fluid 



on an 5^ described by ( 1.1 ) we can readily obtain the tree level retarded Green's function 



of (see for example |43j ) . 

To perform the analysis, it is in fact more convenient to use the (r, p)-coordinate 
system introduced in (2.5). To study different timelike hypersurfaces we fix p = 1 and 
tune a from the horizon at to the worldline at 1/2. We must then study for what 
values of (complex) uj the purely outgoing solutions in (3.10) vanish at the p = 1 



hypersurface. It is relatively straight forward to compute the corrections to the dispersion 
relation perturbatively in a. For instance, to linear order in a we findj3 

l + "(5 + 34) . (4.2) 



Expression (4.2) is only reliable for / 



2 

2/3a. 



4.2 Numerical Results 



r^= 1.0000 X 10~V 



r^= 0.1800 <■ 



r^= 0.3600 <■ 




= 0.5400 1 



r,= 0.7200/ 



r, = 0.9000 /■ 
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Fig. 4: Flow of frequency spectrum iool vs / as we move away from rc ~ toward the 
cosmological horizon. 

Since we impose Dirichlet boundary conditions at p = 1, a parametrizes the location 
of our cutoff surface rc with respect to the cosmological horizon. The relation is given by 

t-Tc l-X 

a = — = . (4.3) 

2^ 2 ^ ' 



^^It is amusing to note that such a correction could be obtained by adding a suitable forcing term to 
our incompressible Navier-Stokes equation [43]. 



13 



As we move Tc away from the cosmological horizon, we expect to deviate from our 
quadratic dispersion relation (3.12). Generically when searching for zeros of Sgf^*{p = 1) 



in (3.10) for arbitrary but fixed rc and /, one runs the risk of finding any one of a tower 



of such zeros (see (4.1) for example). 



r,= 0.9646/ 



.3 60 



r.= 0.9677 r 



r,.= 0.9708; 



I 

rr= 0.9740/ 



/ 



0.9834/ 



0.9928/ 
^ . — J — , — , 



I 

0.9771 



/ 

0.9803/ 




r,= 0.9960/ 



0.9991/ 



"§ 60 - 



Fig. 5: Flow of frequency spectrum iojl vs I as we approach Vc = i- Points obeying the 
fluid dispersion relation are represented in blue. 

In order to make our analysis clear, we perform a search for the lowest lying zeros at a 
given Tc and I and do not present the rest of the tower. In the plots in figures |4] and [5] we 
have not searched for the I = 1 mode as this is where u = and therefore our solutions 
(f"^^ and (/j^" degenerate. We cover this case in appendix [Bj In what follows we will only 
examine the case of pure (negative) imaginary u given that both the quasinormal modes 
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(4.1) and the fluid modes (3.12) obey this property. It would be interesting to extend 
the analysis to general oj in the lower complex plane. 

We now describe how the spectrum behaves as we approach rc = L As a reference, 
we also present the results for rc close to the pole Tc = in figure |4] where the linear 
dispersion relation is apparent. As r^. is increased, we start to observe a deviation from 
the linear behavior and the poles start to cluster into staircase-like behavior. For rc close 
to the horizon, starting from the plot at the top left corner of figure [5} we note that 
there are (at least) three distinct sets of modes separated by large gaps. The / = 2 mode 
lies on the fiuid dispersion relation line (meaning that it satisfies the dispersion relation 
given by (3.12) with quartic corrections as in (4.2)), whereas the rest do not. As we 
move Tc closer to the horizon, we see that the non-fluid modes get pushed higher and in 
the fourth plot, the / = 3 mode is plucked down to the fluid line. This happens once 
again for the I = 4 mode near rc = 0.9865^ whereas the non-fluid modes keep getting 
pulled higher. The reason for these jumps is that the lowest lying zeroes of 5g°^^ are 
modified discontinuously as a function of rc as is visually depicted in figure |8] of appendix 
[C| Finally, we find that arbitrarily close to the horizon the dispersion relation lies entirely 



on the fluid dispersion relation (3.12) computed analytically. 




Fig. 6: Flow of fiuid mode frequencies with fit lines. The gray line is given by iool = 
+ 1) — 2) while the orange line is given by ioJ^ = v{ot) + 1) — 2) with v{a) given 
by (4.2). Note that the orange line fits the data better away from rc = £ and the gray 



and orange lines coincide at the horizon, as expected. 



As a check of our analysis in the previous section, in figure [6] we show that the 
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correction at 0{a) of the fluid viscosity gives the correct behavior for the fluid mode 
frequencies for small but nonzero a. Notice that the orange line in figure |6] fits the data 
better than the gray line for Vc parametrically displaced from the cosmological horizon. 
The orange line is precisely the dispersion relation corrected to 0{a) in our analysis of 
linear perturbations of the background metric of the static patch. The gray line does not 
include 0{a) corrections. 

It is interesting to note that all observed spectra of iuj£ are monotonically increasing 
functions of /. This is a feature that we may expect by continuity away from the modes 
analytically computed at rc = i, but holds as we push rc a finite amount from the surface, 
even when jumps occur. 



5 Incompressible Fluids on Spacelike Slices? 

So far we have discussed several aspects of the static patch, which is the region accessible 
to a single observer. We would now like to briefly discuss some aspects of the future 
diamond. After all, future infinity is clearly a viable candidate location for the non- 
perturbative definition of an asymptotically de Sitter universe. Ordinarily, we would not 
associate the dual theory on with the static patch observer. On the other hand, 
if we impose boundary conditions where there is no incoming flux from the Northern 
diamond, such that all data reaching is coming from a single static patch one might 



conceive of such a relation In this section we will make some simple mathematical 
observations about the behavior of metric deformations on a spacelike slice just outside 
the cosmological horizon. 



5.1 Linearized Analysis 

Our aim is to solve the linearized equations in the future diamond, imposing Dirichlet 
boundary conditions on a constant r surface arbitrarily close to the cosmological horizon. 
In order to choose an appropriate near horizon coordinate system we begin with the 



future diamond, described by (2.2) with r G [£, oo]. As before, we introduce the following 
coordinate transformation: 



t _ 1 , 1 
I ~ 2^^~ 2 



f - ^log(^-p(l-ap)"^) , r = e{l-2ap). (5.1) 



In what follows we will drop the tildes. The metric becomes: 

- ^ +p'^)dT'^ + 2dT dp + {1-2 pa fdnl . (5.2) 



£2 



a 



Taking a — t- with a > 0, constant p slices for p < now correspond to spacelike slices 
just above the cosmological horizon. This slice receives data from the future horizons of 



^^Such 'holographic projections' of the static patch observer were also considered in the rotating Nariai 
geometry 44 , _45i 36). In that case a near cosmological horizon limit allowed for an isolated space-time 
whose (spacelike) boundary is of the same type as the spacelike slice we are discussing here. 
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both the Southern and Northern patches. If we are to isolate one of the observers, say 
the Southern observer, we must impose that the incoming flux from the Northern static 
patch vanishes, as shown in figure [7| 




Fig. 7: Our boundary conditions for the linearized modes are such that the induced 
metric on the p = —1 slice is unchanged and there is no incoming flux from the Northern 
static patch. 

The vector mode with vanishing flux from the Northern patch is generated by a 
master fleld $^ = e-2^"^*(^.p)/^(^^, with: 



where 



-pY 



iFi [ai,bv,ci;ap{-l + 2ap)-^] (1 - ap)-^°'^(l - 2ap) 



\2iauj 



ai = —I — 2iau) , 6i = 1 + / — 2iauj , ci = 1 — 2iauj . 
The linearized metric deformation is given by: 

(1 - 2ap) 



^9: 



<5/ 



2Vi X e-^'^^(-y9)*"'^+i(l - ap)-''^^+^ 
1 — ap 



ip 



-2aVi X e 



-P 



1 



1 



lUJ 



2a 



dp ^: 



s 

V 1 



(5.3) 
(5.4) 

(5.5) 



{l-2ap)\ {l-2ap) 



2p{l - ap) J 



2a 



(5.6) 



Demanding that the above vector modes vanish at the spacelike p = —1 hypersurface 
in the limit where a — )• gives the discrete relation: 



,f = +i{l{l + l)-2) , 1 = 1,2,... 



(5.7) 



This is precisely the same dispersion relation that was found earlier in (|3.12|) for the 
Lorentzian hypersurface but with Vgf = — 1 
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To linear order in a we flnd: 



i^sf = -1 + 1 (5 + 34) 



(5. 



^^We note that a time-reversal r — r transformation leads to v' -> — and thus introduces a sign 
change to the viscosity term u^f in the non-linear Navier-Stokes equation Thus, one can interpret a 

negative viscosity fluid as a time-reversed version of a positive viscosity fluid. We thank R. Loganayagam 
for pointing this out to us. 
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It is not hard to show perturbatively in a that i^s/(q) = —v{—ol). Upon evaluating the 
hnearized metric modes on the 'spacelike fluid' frequencies (5.7) we find that the vector 
modes are regular as we approach the p — )• horizon. 



5.2 Non-linear Analysis 

We would now like to perform a non-linear analysis of the metric deformations in an a 
expansion, again in the context of the spacelike slices. The cutoff surface will now be at 
p = —\. Non- linear deformations analogous to those presented for timelike slicings are: 

IT = (5.9) 

+ p^dr"^ + 2dTdp + d^l + (1 + p) [v'^dr'^ - 2v^dTdx'] + 2pPdT'^ 



a 



(4/j + 2P) dnl - (1 + p)viVjdx'dx^ 



+ - 1) (V^Vi + RjiV^) drdx' + 2vidpdx' - {v"^ + 2P) drdp + 2{l + p) (f)["^dTdx 
+ (^Ap^dnl + 2glfUx'dp + g^fUx'dx^^ + .... 

On the spacelike slice at p = —1, the internal geometry is conformally equal to 

„2 n 



dsia = l^- + I j dr-" + {1 + 2ay dni , (5.10) 
with a conformal factor 

l-2aP + 0{a'^) . (5.11) 

Similarly to the timelike case, for this metric to solve the Einstein equation through 
0{a^), {vi,P) are required to solve the incompressible Navier-Stokes equation: 

drv' + VjV^g^v' + %2P - Vsf i^hv' + Py) = (5.12) 

with 

vsf = -1 . (5.13) 

Note again that the viscosity v^f has changed sign from the fluid on the timelike slices. By 
integrating the QrT = condition at 0{a^) over the sphere, we again find the constraint 
that 

dr( I v'^d^'A = -2dr ( f Pdn2] ■ (5.14) 



Hence the structure of the Einstein equation with positive cosmological constant on the 
timelike slice near the cosmological horizon with the specified boundary conditions is 
closely mimicked in this spacelike context. 
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5.3 Pushing the Spacehke Slice to X+ 

We now wish to push the spacehke shce all the way to X"*" and study the constraints on u. 
We impose fast-falling Dirichlet boundary conditions at and no incoming flux from 
the Northern patch. 

Reverting to static patch (r, t)-coordinates, the solutions ^v{r,t) = e~^'^*'(pv{r) anal- 
ogous to (3.6) and (3.7) near were computed in |14j and are given by: 

-iu}£/2 



r 

^2 



1 



2-Fi ai;6i;ci; 



-iuil/2 



2F1 a2;62;c2; 



with 



ai 



02 



{l + l 
{l-l 



lUJl 



lUJt 



bi 
b2 



{—I — iojt) , c\ 
{2 + 1- iul) 



C2 



(5.15) 
(5.16) 

(5.17) 
(5.18) 

00 we find 99" ~ 1 and ip"^ ~ l/r. Note that 
(c2 —02 — 62) = i^i- In order to eliminate deformations of the conformal metric at X"*" 
we switch off the slow falling mode ip~ . 

Our task is to eliminate the incoming Northern flux upon turning on (/?^. It is not 
hard to see that this will require (c2 — 02 — 62) = to be an integer. We must further 
ensure that metric fluctuations are analytic for iui G Z"*". To achieve this, we must 
eliminate the logarithmic term in the hypergeometric identity (D.l). This implies that 
either 

62 = -n2, ni = 0,1,2,... (5.19) 



As we approach X+ in the limit r 



0-2 



-ni 



or 



It turns out that of the two possibilities, only the second one is sufficient to eliminate the 
Northern incoming flux. In the flrst case, we have to impose a further inequality ni > I, 
whose origin is explained in appendix |Dj Hence, deflning rii = I + hi, hi = 0, 1, 2, . . . 
and imposing no further condition on the integer ?t-2, we obtain the following conditions: 



-i{2hi + 1 + 1) 
-i(2n2 + 2 + 1) 



ni 
n2 



0,1,2,, 
0,1,2,, 



(5.20) 
(5.21) 



which combines into one single tower of modes 

= -i(n + / + 1) , 



n 



0,1,2,., 



(5.22) 



Curiously, and perhaps interestingly, this is exactly the same set of modes as the quasi- 
normal mode spectrum (4.1) in the Southern patch. 
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5.4 Topological Black Holes in AdS4 



In fact, the above calculation is mathematically equivalent to the computation of quasi- 
normal modes for the massless topological black hole in AdS4 [HI HHl [l9] (see also [50] ) 
upon continuing ky — )■ —ky. This is due to the fact that the metric of the massless 
topological AdS4: 

is related to the static patch metric by an analytic continuation. The two-dimensional 
space: d'H2 = {d^,"^ + sinh^^d^^) is the standard metric on the hyperbolic two-manifold. 



The analytic continuation from the static patch metric (2.2) to the above metric is given 
by: 

£ ilAdS , t it , r ^ ir , 9 ^ , (f) ^ (j) . (5-24) 

An observer in the massless topological AdS4 geometry observes a Hawking temperature 
given by T = l/27r^^rf5. If one considers compact quotient of 'H2-, there is a finite 
entropy proportional to iJlAdsl^p)^ associated with the horizon at f = iAdS- We also 
expect such mathematical similarities to hold between the boundary correlators near 
the boundary of topological AdS4 black holes and those near X+ (using the boundary 
conditions discussed above). It is also interesting to note that at the non-linear level one 



can add negative energy to (5.23) and create spherically symmetric asymptotically AdS4 
topological black holes. This occurs up to a critical negative mass, for which one finds a 
solution interpolating between AdS4 near the boundary and AdS2 x 07] . Similarly, 
adding sufficient mass to the static patch leads to the Nariai solution which interpolates 
between dS4 near X"*" and dS2 x 5^. We hope to further uncover this map in future work, 
as it may provide insight into the nature of the static patch observer. 
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A Scalar Perturbations 

Gravitational perturbations in static dS4 consist of a scalar and vector mode of the SO (3) 
of the two-sphere. There are no tensor modes in four dimensions. Scalar harmonic 
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perturbations can be reduced to a single Ishibashi-Kodama master field \37\ l38] , which 
obeys the same effective equation as that of the vector perturbations except that the 
angular number I begins at Z = (instead of Z = 1 in the vector case). 

An incompressible fluid requires a divergenceless velocity field Vi . The scalar harmonic 
allows only the possibility Vi (x Si = — Vj5. S is the scalar harmonic on the sphere which 
satisfies: 

(V|2+A;|)5 = 0, kl = l{l + l), / = 0,1,2,... (A.l) 
Imposing incompressibility leads to 

%2Vi oc VI2S = -klS (A.2) 

which vanishes for ks = 0, i.e. the spherically symmetric mode / = 0. In this case Vi = 
and we are left with a trivial fluid. 

Thus, in the case of an incompressible fluid, the scalar mode only consists of trivial 
fluids. It would be interesting to investigate the case of a compressible fluid which allows 
for sound modes 



B The / = 1 Vector Perturbation 

Gravitational vector perturbations with I = 1 differ from I > 2. We follow the discussion 
in [371 [38]. In addition to equation 3.3, spherical vector harmonics satisfy the following 
equation: 

{Vl2 + k^-3)V^J=0 , Vij = -^[DiVj + DJV^] , (B.l) 

where Di = (V52)^. For /cy — 3 < it can be shown that Vij vanishes and therefore, V* 
must be a Killing vector on the sphere. In this case, we parametrize the perturbations 
as: 

Sgai = rfaVi, (B.2) 

where = {t, r} and = {9, (j)}. Given that Vij = implies Sgij = 0, we can no longer 
fix the gauge freedom by imposing 6gij = 0. Instead, we will fix the gauge = 0. Prom 
[38] . using the only gauge invariant object Fab- 

r-'Fab = Da - Db , Da ^ (V,(2))^ (B.3) 



the Einstein equation imply that Fab is: 



3J 

Fab = -^eab (B.4) 



with J constant. In the case of a perturbation about a spherically symmetric black hole, 
this corresponds to giving it a small amount of angular momentum [38J . 



^''in [20] , it was shown for planar horizons that the speed of sound for the scalar sound mode goes to 
infinity and effectively decouples from the non-relativistic fluid sector. 
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We can study the r dependence of the metric components. Given that the connection 
term will drop out due to anti-symmetrization, we can replace the covariant derivatives 
with ordinary derivatives and find: 



The above integrates to 



3J 



Mr) = J 



(B.5) 



(B.6) 



where we have set Dirichlet boundary conditions on dgu at r = Tc- Hence, the perturbed 
metric reads 



ds = dsQ + J 
with Vi a Killing vector on S'^. 



1 



„2n 



Vi{e,(l))dtdx' 



(B.7) 



C Mind the Gap 

Looking at the last two plots in figure [5j we see that as we move the cutoff surface away 
from the cosmological horizon, the zeroes for I = 9 and I = 10 modes jump discontinu- 
ously. We would like to describe how this behavior arises. Recall that our goal was to 
find the lowest lying zero of Sgf^^{p = 1) as a function of oj on the negative imaginary axis 
as we move re- Notice in figure [8] that as we move Tc to smaller values, the lowest lying 
zero disappears and the new lowest lying zero is at a finite distance away in frequency 
space. 

D Hypergeometric Gymnastics 

For (c — a — b) positive integer the following relation holds |51j : 

/ , ^ r(c - a - 6)r(c) ^ ^ {a)nib)n .n 

T{c-b)r{c-a) ^ [1 + a + b - c)nn\ 

+ - ^>°-°"f(l^g i;:+'y-:-°t)! '^" - - - • (^■'> 

where {a)n = r(a + n)/T{a) are the Pochhammer symbols and 

kn = Tpin + 1) + ip{n + l + c-a-b)- Tp{n + c - a) - il){n + c-b) (D.2) 

with -0(2;) = d\ogT{z) / dz. 

We treat the case covered by equation (D.l ) as it is relevant to the text. To eliminate 



the log terms, we require l/r(a)r(6) = 0. In the case of a = — ni, ni = 0,1,2,..., 
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=0.999/ r, =0.998/ 




Fig. 8: Flow of Sgf^^{p = 1) for / = 10 as we move Tc away from the cosmological horizon. 
We wish to find frequencies that make this function vanishes. As we vary rc the zero 
jumps. 



if c — 6 > (which in the spacehke case becomes ni > I), then the whole second sum 
vanishes and we get [5.2]: 

2F1 (a; b; c; z) = (a; b;l-c;l-z) (D.3) 

i (c — oji (c — a) 

which goes to a constant as r — )• 1 (here z = £^/r^) and translates into the modes with 



no incoming flux from the Northern patch. So, these are the modes we want,^^ 

In the other regime c — b < 0, the first sum vanishes completely due to the gamma 
function's poles. Naively, we would think that the whole expression is zero, however, the 
second term contains a divergent term ip{n + c — b) in the sum that will cancel out the 
divergence in l/r(a) for n < — (c — 6). then we actually have 



^^In the case of 6 = —n2, n2 — 0, 1, 2, . . ., the analogous inequality in the first regime is c — a > 0. 
However, our parameters imply that c — a > is always true. Thus we are always in this first regime 
which gives the modes that we want and there is no further restriction on 722. 
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leading to[52]: 

2F1 (a; b; c; z) ~ r(c)r(-Q + l) _ ^ , _ ^. ^ _ ^ + ^ _ ^) , (D,5) 

(c — a — oj!l [b) 

This implies that the 99^ tends to (r^ — xy^^"^ (as r — )• 1) which is an incoming wave 
from the Northern patch and should be excluded. 
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